Nonequilibrium critical dynamics in inhomogeneous systems 
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We study nonequilibrium dynamical properties of inhomogeneous systems, in particular at a free 
surface or at a defect plane. Thereby we consider nonconserved (model-A) dynamics of a system 
which is prepared in the high-temperature phase and quenched into the critical point. Using Monte 
Carlo simulations we measure single spin relaxation and autocorrelations, as well as manifold auto- 
correlations and persistence. We show that, depending on the decay of critical static correlations, 
the short time dynamics can be of two kinds. For slow decay of local correlations the usual domain 
growth process takes place with non-stationary and algebraic dynamical correlations. If, however, 
the local correlations decay sufficiently rapidly we have the so called cluster dissolution scenario, 
in which case short time dynamical correlations are stationary and have a universal stretched ex- 
ponential form. This latter phenomenon takes place in the surface of the three-dimensional Ising 
model and should be observable in real ferromagnets. 



I. INTRODUCTION 

To characterize the dynamical universality class of a 
critical system at equilibrium it is enough to provide 
the value of the dynamical exponent, z, which generally 
depends on the local dynamics, conservation laws and 
symmetries J* In a nonequilibrium system, which is pre- 
pared by quenching it from the high-temperature initial 
state to the critical point, new nonequilibrium critical 
exponents have to be defined. The reason of this is the 
broken time translation invariance due to a discontinuity 
at the time horizon ("time-surface"). Since at the critical 
point there is no characteristic time-scale the nonequilib- 
rium preparation has a long-lasting effect, which has two 
important consequences. First, the dynamical correla- 
tions, such as the single spin autocorrelation function, 
C(t,s), are generally non-stationary. They depend on 
both the preparation or waiting time, s, and the obser- 
vation time, t. The second effect is the appearance of a 
new nonequilibrium critical exponent, A, which is used 
to describe the long-time asymptotic, t 3> s, of the auto- 
correlation function as C(t, s) ~ t~ x / z X 

Another, and related process of nonequilibrium dy- 
namics is the relaxation of the magnetization, if in the 
initial state there is a small, non- vanishing value of mi. 
For short times the magnetization has a power-law de- 
pendence: m(t) ~ mil , with the initial slip exponent 8 
which satisfies the scaling relation: A = d — dzS^ Here 
d is the spatial dimension of the system. In mean-field 
theory 8 = 0, whereas in real system generally 8 > 0. 
Thus due to fluctuations the order is increasing in the 
early time regime. This phenomenon is related to the 
fact that in the initial state there are no long-range cor- 
relations, thus the systems is mean-field like. Since the 
actual critical temperature of the system, T c , is lower 
than the critical temperature of the mean-field model, 
T c (mf) > T c , in the early time-steps there is an effec- 
tive coarsening process, during which the magnetization 
is increasing. The nonequilibrium magnetization with 



rrii < 1, however, can not exceed the value obtained 
with m t = 1, when from the ordered initial state the 
magnetization decreases as m(t) ~ t~ x l z , where x is the 
anomalous dimension of the magnetization. This criti- 
cal parameter is defined trough the asymptotic decay of 
the equilibrium, equal-time correlation function at the 
critical point: (ai(t)<jj(t)) ~ \i — j\~ 2x and satisfies the 
scaling relation: x = (3/v, with [3 and v being the usual 
static magnetization and correlation length exponents, 
respectively. From the previous reasoning we obtain for 
the border of the short-time regime, t$, as ti ~ m i z ' x, j 
where Xi — 8z + x is the anomalous dimension of the ini- 
tial magnetization. Note that the nonequilibrium time- 
scale, U, diverges as m; — > 0. 

Besides single spin autocorrelations one often consid- 
ers the dynamics of extended objects, such as the order- 
parameter of the whole sample (global autocorrelations) 
or that of a manifold 5 with dimension d! < d. In the 
asymptotic regime t 3> s the manifold autocorrelation 
function displays an algebraic decay 5 -: G(t) ~ i~ A l z with 
A' = A — d' . Still another dynamical quantity is repre- 
sented by the persistence, P pr (t), which measures the 
probability that the order-parameter associated either to 
a single spin, or to a manifold, or to the whole sample 
has not changed sign within time t& Among others, per- 
sistence is related to the properties of the given autocor- 
relation function and for global and manifold persistence 
it has often a power-law dependence: P pr (t) ~ t _e " r , 
where Q pr is a new nonequilibrium exponent X Only for 
global persistence and for a Markovian process can Q pr 
be expressed by other known exponents^ 

The dynamical critical behavior, as outlined above, 
refers to the bulk of a homogeneous system in the thermo- 
dynamic limit. Real materials, however, are bounded by 
surfaces and many nonequilibrium processes (thcrmaliza- 
tion, transfer of heat, etc.) are very intensive at the sur- 
face. At a free boundary layer the order is weaker than in 
the bulk due to missing bonds and correlations between 
two surface sites, (cr^ (t)aj 1 (i)) ~ \ii — ji\~ 2xi in- 
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volve a new scaling dimension, x\ > ^^iiiSiiiia This 
satisfies the scaling relation: x\ — @i/v, with f3% being 
the surface magnetization exponent. Note that there is 
an analogy between spatial and temporal surfaces, where 
translational invariance in space and in time, respec- 
tively, is broken. Similarly the scaling dimensions, x\ 
and Xi, respectively, play analogous roles. According to 
field-theoretical investigations^ no new dynamical criti- 
cal exponent should be introduced at the surface: z and 
Xi remain unchanged whereas Q\ and A l5 as defined on 
surface spins, can again be expressed by the known ex- 
ponents, including xi^ii 

We note that besides a free surface there exist other 
types of inhomogeneities, such as a localized or an ex- 
tended defect plane, near which scaling of the local mag- 
netization is different from that in the bulk. Correlations 
between two sites at the defect plane, ii,ji are asymp- 
totically given at the critical point by (c^, (t)tTj, (t)) ~ 
\ii — ji\~ 2x ' , which involves a new static (local) expo- 
nent Xi 7^ xi& In an inhomogeneous system it can be 
found in many cases that xi > Xi, thus the spatial sur- 
face or interface is more disordered than the "temporal 
surface" . This happens for example at the ordinary sur- 
face transition of the three-dimensional Ising model, with 
x\ = 1.26 and x\ — 0.74, and should be realized in real 
magnets, too. Nonequilibrium dynamical critical behav- 
ior of systems with xi > Xi is expected to be different 
from that with xi < a;,. In the latter case traditional 
domain growth (DG) scenario should hold. On the other 
hand in the former case critical relaxation is more effec- 
tive than the nonequilibrium growth, therefore no new 
domains are created in the surface or interface region. 
The dynamical process is dominated by rare clusters with 
correlated sites, and the nonequilibrium dynamics is re- 
alized through cluster dissolution (CD). The CD process 
can be seen also for the surface manifold autocorrelation 
function, provided d! — xi — x < 

Our aim in the present paper is to study in detail the 
nonequilibrium dynamical critical behavior in inhomoge- 
neous systems, in particular at a free surface or at de- 
fect planes. Here we restrict ourselves to nonconserved 
(model-A) dynamicsi We use scaling theory and perform 
extensive Monte Carlo simulations. We clarify the prop- 
erties of the nonequilibrium dynamics in the DG and CD 
regimes. A short account of our results obtained at a free 
surface has been announced in a Letter m& 

Our paper is organized as follows. In Sec. 2 we 
start with the presentation of a phenomenological picture 
based on scaling theory. Sec. 3 contains our numerical re- 
sults on local magnetization relaxation and autocorrela- 
tions i) at free surfaces in the two- and three-dimensional 
Ising models, ii) at an extended surface defect in the two- 
dimensional Ising model (Hilhorst-van Lccuwcn (HvL) 
model) ^* and iii) at an internal defect line in the two- 
dimensional Ising model (Bariev model) m& We also con- 
sider iv) surface manifold autocorrelations and persis- 
tence. In Sec. 4 we close our paper by a discussion. 



II. SCALING THEORY OF NONEQUILIBRIUM 
DYNAMICS IN INHOMOGENEOUS SYSTEMS 

We consider a d-dimensional system having a free sur- 
face or a defect plane, the distance of it being denoted 
by y. (In the following we will refer to that plane as 
a surface and the local scaling dimension will be called 
xi, so that xi — x\ at a free surface.) From the high- 
temperature phase, T > T c , we quench the system to its 
critical temperature and study different physical quanti- 
ties (magnetization, single spin autocorrelation, manifold 
autocorrelation) as a function of the waiting time, s, and 
the observation time, t, both measured from the time of 
the quench. We start our analysis with the relaxation of 
the magnetization. 



A. Relaxation of the magnetization 

In this case the system is prepared with a small ini- 
tial magnetization, m, > 0, and the local magnetization 
m = {<J y (t)) = m(y, t, m») is measured at the critical tem- 
perature. Here a y {t) is the operator of the local magne- 
tization. Under a scaling transformation, I —> l/b, when 
length, /, is rescaled by a factor, b > 1, the magnetization 
behaves as: 

m{y,t,mi) = b- x m(y/b,t/b z ,m l b^) . (1) 

Here Xi is the new, nonequilibrium scaling dimension of 
the initial magnetization, as introduced in the previous 
section. Taking in Eq.QJ b — t 1 ^ , we arrive to: 

m(y, t, rm) = r^'miy/i 1 '", m^l z ) , (2) 

in which the scaling function, m(r, fx), behaves differently 
in the bulk, r — > oo, and at the surface, r — * 0, respec- 
tively. 



1. Bulk behavior 

In the bulk of the system we have by definition 
linv^oo m(r, fx) = rhb(fi), and the bulk scaling function 
mb(/i) has a cross-over: 

For short times, t < ti ~ m i z ^ Xz , the magnetization in- 
creases as 7Ti (i) ~ t e , with 9 — (xi — x)/z. For longer 
times, however, the magnetization follows the equilib- 
rium decay, m(t) ~ t~ x l z . We note that the short-time 
behavior of the magnetization is the result of two differ- 
ent processes. Due to nonequilibrium domain growth the 
magnetization increases as t Xi / z , whereas critical relax- 
ation acts to reduce its value by t~ x / z . 
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2. Surface behavior 

In order to understand the concept of surface critical 
behavior we start with the equilibrium magnetization, 
m e q(y,t,6), where S = (T c — T)/T c is the reduced critical 
temperature. The scaling transformation now reads as: 

m cq (y, t, 6) = b- x m cq (y/b, t/b z , 5b 1 '") . (4) 

We consider first the static behavior with t = 0, set the 
length-scale to b = 5~ v and obtain m eq (y,t = 0, S) = 
5 x "rh eq (jj5 u ). Here we notice the temperature depen- 
dence of the bulk magnetization: m^ q ~ 8@ , with a 
critical exponent [3 = xv and with the scaling function: 
linif^oo m cq (f ) = const. By definition the magnetization 
at the surface behaves as: m\ q ~ 5" 1 with /3; = x\V. This 
is compatible with the above scaling result provided the 
scaling function for a small argument (short distance) be- 
haves as: Iiruf— >o m oq (r) ~ f Xl ~ x . Now considering the 
equilibrium dynamical behavior at the critical point we 
set in Eq.QJ 5 = and b — t 1 ^. Repeating the previous 
reasoning we obtain for the magnetization in the surface 
region m cq (y,t) = t~ x ' z fh eq (r) ~ t~ x / z r Xl ~ x , r -C 1, 
which can be obtained from the short-distance expan- 
sion in field-theoretical calculations^ Note that starting 
from the bulk order is decreasing towards the surface for 
xi > x. In nonequilibrium relaxation two processes have 
to be taken into account close to a surface: the magneti- 
zation i) increases due to nonequilibrium domain growth 
and ii) decreases due to surface effects. Of course, the 
first process is of no importance when starting from a 
fully ordered statei^LSSi The relative strength of the two 
processes are given by the scaling dimensions x% and xi, 
respectively. Nonequilibrium dynamics is markedly dif- 
ferent for xi > xi, which is the domain growth regime, 
and for Xi < Xi, which is the regime of cluster dissolution. 

a. Domain growth: Xi > x\ In this case bulk 
nonequilibrium order penetrates into the surface region, 
in which, as time goes on, new domains are formed. 
Therefore short distance expansion of the magnetization 
stays valid, yielding for small r the following expression 
for the scaling function m(r, /i) given in Eq.Q: 

lim rh(r, /x) = r Xl ~ x fhi(u), x\ > x\ . (5) 

Here to/(/i) has the same limiting /z-dependence as fhb(p) 
in Eq.Q. Thus we have in the short-time limit, t < t j ~ 
mJ z ' Xi , mi(t) ~ t 6 ', with 0i = 0—{xi~x)/z = (x l -xi)/z. 
Thus the initial slip exponent at the surface in the DG 
regime can be expressed by known exponents^ 

b. Cluster dissolution: Xi < x\ In this case, for short 
times nonequilibrium order can not penetrate into the 
surface region and therefore no new domains are formed 
there. In contrary the initially existing ordered clusters 
are diluted and the initial relaxation has a fast, non- 
algebraic dependence. Evidently the short distance ex- 
pansion is not valid here. Analyzing the CD process we 
have determined in the Appendix the small r behavior 



of the scaling function, m(r, fi). It is given in a stretched 
exponential form: 

lim rh(r, /i) ~ cxp(— Cr K ), Xi < xi , (6) 

r— >0 

with a universal power: 

_ (gj - x{)d 
K ' d-l ' (7) 

If k/z > 1, as argued in the Appendix, the functional 
form in Eq.© is pure exponential. This time dependence 
is valid for t < t[ ~ | In77ii| 2 / K , for larger times one goes 
over to the DG regime with the decay in Eq.©. 

B. Single spin autocorrelation function 

The nonequilibrium autocorrelation function is defined 
by C(y,t,s) — (o- y (t)a y (s)) and in the prepared state 
there is no initial magnetization. From the scaling trans- 
formation of the autocorrelation function: 

C(y,t, S ) = b- 2x C(y/b,t/b z ,s/b z ) (8) 

we obtain with b = t 1 / 2 : 

C(y,t,s)=t- 2x / z C(y/t 1 / z ,s/t) . (9) 

In the following we analyze the scaling function, C (r, r) , 
t < 1, in the different limits. 

1. Bulk behavior 

The scaling function in the bulk, linv^oo C(r, r) = 
Cb(r), has the asymptotic behavior: 

lim (%(t) ~ T W-«i-*)/* . (10) 

So that for t > s we have Cb(t) ~ t~ x l z , with the bulk 
autocorrelation exponent: A = d — x^ + x. 

2. Surface behavior 

As for the magnetization relaxation we have to dis- 
criminate between two regimes: the DG and the CD dy- 
namics, for Xi > xi and Xi < xi, respectively. 

a. Domain growth: Xi > xi In this case domain 
growth takes place in the surface region, too. There- 
fore for small r we can separate the r-dependence from 
the scaling function in Eq.JSJl, leading to: 

lim C(r, t) = r^-^C^r), a* > s, . (11) 

r— >Q 

The small r dependence of the scaling function C'i (t) is 
the same as in Eq. l(TU)l for the bulk quantity: Ci(t) ~ 
T (d-Xi-x)/z _ c onse q uen tly i n the limit i>swe have for 
the surface autocorrelation function Cj(t) ~ t~ Xl l z , with 
the surface autocorrelation exponent: A; = A+2(x;— x) = 
d — Xi — x + 2x/>i& 
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b. Cluster dissolution: Xi < xi In this case typically 
no new domains are formed in the surface region. The 
autocorrelation function of the system is therefore domi- 
nated by such rare regions in which at preparation there 
is an ordered cluster. The probability of the existence of 
such a region is very small as it decreases exponentially 
with its volume. The relaxation time, however, associ- 
ated to such a cluster is very large, it is the exponential 
function of the typical size of an interface in this clus- 
ter. Combining these two effects we arrive to a stretched 
exponential dependence, as described in the Appendix. 
Since the starting cluster structure does not depend on 
the waiting time the surface autocorrelation function for 
a small t — s <C s is stationary. In this case both r — > 
and 1 — r — > 0, so that the appropriate scaling combi- 
nation is p = (1 — r) 1 / 2 jr. In terms of this the surface 
autocorrelation function in the CD regime is 

<7,(t)~exp(-C"(t -«)"/*) (12) 

for k/z < 1, with the exponent k given in Eq. J7J), 
whereas for k/z > 1 the functional form in Eq. I|12l) is 
pure exponential. 



C. Manifold autocorrelations 

Here we consider a d! dimensional manifold, M, which 
is located at a distance y from the surface. The mani- 
fold order-parameter is the sum of the individual opera- 
tors, S v = J^ieM °V(0j an d the manifold autocorrelation 
function is denned by G(y,t,s) = (S y (t)S y (s)) . This 
quantity, as well as manifold persistence, has been intro- 
duced and studied in the bulk of the system by Majumdar 
and Bray£ First, we consider the equal-time autocorrela- 
tion function, G(y,t,t) = (J2ieM (J y ( S)T, ] eM °v(?))> and 
note that in the disordered initial state at t = s = it 
is zero, since S y = 0. As time goes on correlations be- 
tween sites within the correlated domain of linear size, 
£ ~ t 1 ^, are built, so that for a bulk manifold con- 
tribution of terms with alternating signs will lead to 
Wyi^T^jeM^U)) ~ t~ 2x / z and therefore G b (t,t) = 
lim^oo G(y, t, t) ~ ^ d '- 2x ~>/ z , provided d' > 2x. If, 
however, d' < 2x, since the equal-time autocorrelation 
function can not decrease, Gb{t,t) will approach a finite 
limiting value within a characteristic (microscopic) time. 
In this case G&(i, t) being independent of the waiting 
time the manifold autocorrelation function is stationary, 
Gb(t, s) = Gb(t — s). This type of results are in complete 
agreement with the mean-field and scaling results of Ref. 

The above observations can be cast into a scaling re- 
lation for the manifold autocorrelation function as: 

G(y,t,s) =t( d '- 2x V*G(y/t 1 / z ,s/t) , (13) 

in which the scaling function, G(r, r), has different ana- 
lyticity properties in the different regimes. 



1. Bulk behavior 

The scaling function in the bulk is denoted as 
lim r ._ >0 c G(r, r) = G 6 (t). 

a. N on- stationary regime: d' > 2x Here G&(t) is 
analytic for r > 1 and has the limiting behavior: 

lim G 6 (r) ~ ,-(«*-*«-*)/* . (14) 

So that for t > s we have Gb(t) ~ t~ x / z , with the bulk 
manifold autocorrelation exponent: A' = d — d! — Xi + x. 
Note that for the global autocorrelation function with 
d' = d the decay exponent, X' /z, is just the initial slip 
exponent, 8. 

b. Stationary regime: d' < 2x In this regime the 
scaling function is non-analytic at r = 1 and given for 
1 - t <C 1 as: 

Gb(T)^(l-T + r Y d '- 2 ^/\ (15) 

where tq is related to the microscopic time-scale until 
which Gb{t,t) reaches its limiting value. In the limit, 
1 — t tq , we combine Ea. l|15|l with i|13|) and obtain a 
stationary short-time behavior: 

G b (M) ~ (i ~ s) (d '~ 2x)/z , (16) 
in complete agreement with Ref. 0. 

2. Surface behavior 

Here we consider first the initial value of the manifold 
autocorrelations at the waiting time t — s, thus we have 
r = 1. For equal time operators the short distance expan- 
sion is expected to be the same as for the magnetization 
in Eq.© so that the scaling function behaves as: 

lim G(r, t = 1) - r x '- x . (17) 

r— >0 

Consequently the surface manifold autocorrelation func- 
tion at the waiting time reads as: 

G t (t,t) . (18) 

Thus, as for the bulk manifold, the time-dependence of 
Gi(t,t), and thus the non-equilibrium dynamics has two 
different regimes. For d! > xi + x, when Gi(t, t) is diver- 
gent for large t we are in the standard DG regime. On 
the contrary for d' < xi + x, when Gi(t,t) approaches a 
finite limiting value within a microscopic time-scale the 
dynamics for t > s is of the CD type. 

a. Domain growth: d' > xi + x For t > s the mani- 
fold autocorrelation function has the same short distance 
expansion as the single spin autocorrelation function in 
Ea. (|ll|) . thus the scaling function, G(r,r), for small r 
behaves as: 

lim G(r,r) = r 2< - x '~ x) Gi(T), d' > x x + x . (19) 
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Note that the form of the short distance expansion is 
radically changing for t > s in comparison with that 
for t = s in Ea. (|17|l . The scaling function, Gi(t), has 
the same small r dependence as its bulk counterpart 
in Eq.jg, thus lim r ^ G;(r) - T (d-="i-x) / z _ As a fi _ 
nal result for f S> s the surface manifold autocorrelation 
function in the DG regime behaves as: Giit) ~ t~ x '/ z , 
with the surface manifold autocorrelation exponent: A[ = 
A' + 2 (a;; — x) — d — dl — Xi — x + 2xi. 

b. Cluster dissolution: d' < xi + x In this case the 
correlated sites form isolated clusters and the dominant 
contribution of the manifold autocorrelation function is 
due to large, but rare clusters. The dilution of the clus- 
ters during the relaxation process for t > s is similar to 
that of the single spin autocorrelation function. The only 
difference is that here the creation of sites in the domain 
wall of ordered clusters goes on with time as t Xl ~ x , which 
follows from the difference in the short distance expan- 
sions of G(r, t = 1) in Eq. ljT7|) and that of G(r, r < 1) in 
Eq. (|19|l . As for the single spin autocorrelation for short 
times, t — s <C s, the manifold autocorrelation function 
is stationary. This means that as r — > and 1 — r — > 
the appropriate scaling combination is p = (1 — r) 1 / z /r, 
in terms of which the surface manifold autocorrelation 
function is given by 

G,(t) ~exp(-C"(t-sr'/ z ) (20) 

for < k'/z < 1 with the exponent k' — (xi—x)d/(d—l). 
If k'/z > 1 then the functional form in Ea. (|20|l is pure 
exponential. 

Here we note that for manifold autocorrelations there 
is an analogy between the stationary regime in the bulk 
and the CD regime at the surface. The difference in the 
functional forms in Eas.l|16p and 12Ufl is due to the rel- 
evant inhomogeneity in the latter case: bulk order can 
only in a very limited way penetrate into the surface re- 
gion. Formally, in the bulk k' should be replaced by zero 
and therefore the decay is algebraic, as given in Ea. fllGI . 

D. Persistence of manifolds 

Persistence of manifolds is given by the probability 
P pr {t) that the manifold order-parameter, (S y (t)), has 
not changed sign up to time t. The bulk manifold per- 
sistence has been analyzed by Majumdar and Bray^ and 
here we shortly recapitulate their findings. 



1. Bulk behavior 

In this case the functional form of persistence can be 
formally related to the manifold autocorrelation func- 
tion at the waiting time, which according to Ea. (|13|l is 
given by G(t,t) ~ t~ c , with ( = (2x - d')/z. In the 
actual analysis the theorem by Newell and Rosenblatt 23 
is used in which the long time decay of the persistence 



of a Gaussian stationary process is related to the asymp- 
totic decay of its stationary autocorrelation function. For 
C < the non-stationary correlator after rescaling and in 
logi variable is transformed to a stationary one. From 
this persistence in t is obtained in a power-law form, 
P pr (t) ~ i~ e P'', with a non-trivial exponent, Q pr - For 
< £ the autocorrelation function in Ea. H16fl is station- 
ary and one can apply directly the theorem by Newell and 
Rosenblatt i^ 3 . For < ( < 1 persistence is in a stretched 
exponential form: P pr (t) ~ exp(— at''), whereas for £ > 1 
it is pure exponential: P pr {t) ~ exp(— bt). 

2. Surface behavior 

Here we should distinguish between the DG and the 
CD regimes. In the DG regime having a nonstationary 
autocorrelation function persistence can be analyzed in 
a similar way as for the analogous regime in the bulk. In 
this way we obtain that for d' > xi + x, the surface man- 
ifold persistence is of a power-law form, P pr {t) ~ t~®v T . 
In the CD regime, with d' < x\ + x, the autocorrelation 
function is stationary in the short-time regime, t — s < s, 
but for longer times it becomes non-stationary. Since the 
dominant variation of persistence is related to the sta- 
tionary part of the correlator it is natural to assume that 
persistence has the same type of stretched or pure expo- 
nential time dependence, as the autocorrelation function. 
Thus for < k'/z < 1, P pr (t) ~ exp(-ai K '/ z ), whereas 
for k'/z > 1 we have P pr (t) ~ exp(— bt). 

III. NUMERICAL RESULTS 

We discuss in the following the nonequilibrium dynam- 
ical behavior of various inhomogeneous systems and con- 
front the numerical results with the predictions coming 
from the scaling theory just presented. The models stud- 
ied include semi-infinite Ising models in two and three 
dimensions, the Hilhorst-van Leeuwen model (i.e. two- 
dimensional semi-infinite Ising model with an extended 
surface defect) and Ising models with a ladder defect. 

A. The models 

1. Semi-infinite Ising models 

Static and dynamical properties of critical semi-infinite 
systems have been studied very intensively during the last 
thirty yearsii*i2ii^. The simplest of these models is the 
semi-infinite Ising model defined by the Hamiltonian 

U = - J s ^2 a * a i ~ Jb X! (Jl ' T 3 ( 21 ) 

surface bulk 

where the spin located at site i takes on the values 
cTj = ±1. The first sum in Eq. (|21|l is over nearest neigh- 
bor pairs of surface spins, whereas the second sum is over 
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TABLE I: Static and dynamic critical quantities of the two- 
and the three-dimensional Ising models, x: bulk scaling di- 
mension, x\: surface scaling dimension, Xi: scaling dimension 
of the initial magnetization, z: dynamical scaling exponent. 
OT: ordinary transition, SP: special transition point. 





X 


Xl 


Xi 


z 


d = 2 


1/8 


1/2 


0.53 


2.17 


d = 3, OT 


0.516 


1.26 


0.74 


2.04 


d = 3, SP 


0.516 


0.376 


0.74 


2.04 



nearest neighbor pairs where at least one spin is not lo- 
cated at the surface. J s resp. J& is the strength of the 
surface resp. bulk couplings. 

The surface phase diagram of the semi-infinite Ising 
model is well established. It is especially simple in two 
dimensions where surface and bulk order at the bulk 
critical temperature T c , independently of the strength 
of the surface couplings. In three dimensions, however, 
different phase transitions are observed, depending on 
the value of J s . For weak surface couplings order can 
not be maintained at the surface independently of the 
bulk. Consequently, surface and bulk both order at the 
bulk critical temperature. This transition is called ordi- 
nary transition. For strong surface couplings, however, 
with J s > r sp Jb where r sp w 1.50 for the simple cu- 
bic lattice^SiSi the surface can order independently of 
the bulk at a temperature higher than T c . Reducing the 
temperature the bulk than orders at T c in presence of an 
already ordered surface (extraordinary transition). The 
different critical lines meet at the multicritical special 
transition point r sp . 

In our study of the nonequilibrium critical dynamics 
in Ising models with free surfaces we restricted ourselves 
to the ordinary transition and to the special transition 
point, thereby focusing on the short-time critical dynam- 
ics, the dynamical scaling regime encountered for longer 
times having been studied recently in Ref.l24landl25l Ta- 
ble [I] lists the values of the static and dynamic bulk and 
surface exponents for the cases under investigation. For 
the ordinary transition we chosed J s — Jb, but some sim- 
ulations were also done at other values of J s to check for 
universality. As usual, the system was initially prepared 
in an uncorrelated state and then quenched to the criti- 
cal temperature at time t — 0. The dynamical evolution 
of the system was studied using heat-bath dynamics. In 
two dimensions systems with typically 300 2 spins were 
simulated, whereas in three dimensions a typical sam- 
ple contained 60 3 spins. Free surfaces were considered in 
one direction, whereas periodic boundary conditions were 
applied in the remaining directions. No finite-size effects 
were observed for the large systems under investigation. 



2. Hilhorst-van Leeuwen model 

The Hilhorst-van Leeuwen model is a two-dimensional 
semi-infinite Ising model with an extended surface 
defecti 13 ! 19 We have studied this model on a square lat- 
tice. Whereas the couplings parallel to the surface have 
a constant strength J\ , the strength of the couplings per- 
pendicular to the surface vary as a function of the dis- 
tance y to the surface: 

Mv)~ J 2 (oo) = 4- (22) 

y 

with A = A T c sinh(2 J 2 (oo) /T c ) /4i2S, In the present work 
only the marginal case with u — 1 and J\ = J2(oo) has 
been considered. For A < 1 exact results show that the 
local scaling dimension x\ is a continuous function of A 
with 

Xl = \(l-A). (23) 

This remarkable property allows us to study in a very 
systematic way the nonequilibrium dynamical behavior, 
passing from the domain growth regime to the cluster 
dissolution regime by only changing the value of the pa- 
rameter A. Especially, for A = —0.06 we have x\ = Xi, 
whereas for A — we recover the pure semi-infinite model 
with x\ = i . Of course, the nonequilibrium bulk quan- 
tities Xi and z are left unchanged by the presence of an 
extended surface defect. 



3. Bariev model 

The plane Ising model with a defect line is another 
inhomogeneous system displaying non-universal critical 
behavior. In his work Barie\ii^& analyzed two types 
of defect lines: a chain defect, where a column of per- 
turbed couplings with strength J c h is considered, and a 
ladder defect, where modified couplings of strength J; 
connect spins belonging to two neighboring columns (i.e. 
for Ji = the system is separated into two semi-infinite 
parts). The short-time nonequilibrium critical behavior 
has already been studied in systems with a chain defect^ 
In the present work we focus on Ising systems with a lad- 
der defect. 

Bariev's exact results demonstrate the dependence of 
the local scaling dimension xi on the values of the defect 
coupling. For the ladder defect he obtains 

2 

xi = — ^ arctan 2 (ft, ) (24) 

with 

ki = tanh(^)/tanh(^-). (25) 

J is the strength of the unperturbed interactions (in the 
following we set J — 1), whereas T c is the critical temper- 
ature of the pure two-dimensional Ising model. Eq. I|24|) 
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shows that enhanced (reduced) defect couplings yield a 
lower (higher) local scaling dimension as compared to the 
perfect two-dimensional Ising model, the largest value be- 
ing obtained for the semi- infinite model with J; = 0. We 
therefore have always xi < Xi in the case of a ladder 
defect and accordingly expect to encounter the domain 
growth regime for all values of the defect coupling. 

As for the two-dimensional semi-infinite model and the 
Hilhorst-van Leeuwen model systems with the linear ex- 
tent L = 300 have been studied typically. For the Bariev 
model periodic boundary conditions are considered in all 
directions. 



B. Local magnetization relaxation 

In order to study the relaxation of the magnetization 
the system is prepared in an initial state with a small 
magnetization rrii. In the actual simulations initial mag- 
netizations mi — 0.02 and 0.04 were usually considered, 
for some selected cases other values of rrii were also in- 
vestigated. The data shown in the different figures have 
been obtained after averaging over at least 60000 realiza- 
tions using different random numbers and different initial 
states. 

The relaxation of the local magnetization is shown in 
FiguresQ]-|21for the different models. Recall that depend- 
ing on the values of xi and Xi two different scenarios are 
possible. When cluster dissolution takes place [x\ > Xi) 
scaling theory predicts a stretched exponential decay © 
at early times followed by a power-law decay at later 
times governed by the exponent 9i = (xi — xi)/z < 0. 
When xi < Xi, however, the local magnetization should 
display a power-law increase with the exponent 9i = 
(xi — x{)/z > 0. Our data completely agree with this 
picture. 

Looking at Table [I] we see that for the three- 
dimensional semi-infinite model x\ > Xi at the ordinary 
transition, whereas x\ < Xi at the special transition point 
and in the two-dimensional model. Accordingly, we ex- 
pect the cluster dissolution process to take place only for 
the ordinary transition in the three-dimensional system. 
Figure ^ shows the relaxation of the surface magnetiza- 
tion for the pure semi-infinite systems. We first remark 
that at later times a power-law behavior is observed in 
all cases. The values of the exponent extracted in this 
regime are listed in Table HT1 and compared with the ex- 
pected theoretical value Q\ = (xi — x\)j z. In all cases a 
nice agreement is found between the numerical value and 
the theoretical prediction. At early times the surface 
magnetization at the ordinary transition in three dimen- 
sions displays a nonalgebraic behavior. The exponent k 
obtained from fitting the data with a stretched exponen- 
tial (grey line) is listed in Table ITTT1 It nicely agrees with 
the scaling prediction (0). 

The relaxation behavior of the surface magnetization 
is shown in Figure |21 for the Hilhorst-van Leeuwen model 
with various values of the scaling dimension x\. Clearly, 



TABLE II: Numerically determined values of the exponent 0i 
governing the power-law behavior of the local magnetization 
when starting from a state with a non-vanishing initial mag- 
netization. The theoretical prediction is (9; = (xi— x{)/z. OT: 
ordinary transition, SP: special transition point. 



pure senn- 


infinite models 




Xi 


numerical 


theoretical 


et 2 


1/2 


0.015(3) 


0.014 


d = 3, OT 


1.26 


-0.255(4) 


-0.255 


d= 3, SP 


0.376 


0.179(4) 


0.178 


Hilhorst-van Leeuwen model 


A 




numerical 


theoretical 


0.75 


1/8 


0.188(3) 


0.187 


0.50 


1/4 


0.131(3) 


0.130 


0.25 


3/8 


0.076(4) 


0.072 


-0.25 


5/8 


-0.045(2) 


-0.044 


-0.50 


3/4 


-0.107(4) 


-0.101 


-0.75 


7/8 


-0.157(3) 


-0.159 


-1 


1 


-0.214(5) 


-0.217 


-1.25 


9/8 


-0.29(1) 


-0.274 


-1.50 


5/4 


-0.37(3) 


-0.332 


Bariev model 


Ji 


Xl 


numerical 


theoretical 


0.2 


0.376 


0.072(1) 


0.071 


0.4 


0.278 


0.118(1) 


0.116 


0.6 


0.208 


0.148(1) 


0.148 


0.8 


0.159 


0.171(1) 


0.171 


1.0 


0.125 


0.188(1) 


0.187 


1.2 


0.101 


0.199(1) 


0.198 


1.4 


0.085 


0.208(2) 


0.206 



TABLE III: Numerically determined values of the exponent 
k governing the stretched exponential behavior of the cluster 
dissolution process taking place when xi > Xi. The values 
have been extracted from (I) the relaxation of the local mag- 
netization and (II) from the short-time behavior of the single 
spin autocorrelation. The theoretical prediction is given by 
Eq. 0. OT: ordinary transition. 



pure semi-infinite models 




Xl 


I 


II 


theoretical 


d — 3, OT 


1.26 


0.81(3) 


0.81(2) 


0.78 


Hilhorst- 


van Leeuwen model 


A 


Xl 


I 


II 


theoretical 


-0.25 


5/8 


0.11(4) 


0.25(4) 


0.19 


-0.50 


3/4 


0.45(3) 


0.41(3) 


0.44 


-0.75 


7/8 


0.71(3) 


0.65(3) 


0.69 


-1 


1 


0.97(4) 


0.89(3) 


0.94 


-1.25 


9/8 


1.17(4) 


1.19(3) 


1.19 


-1.50 


5/4 


1.48(4) 


1.54(5) 


1.44 
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FIG. 1: Relaxation of the surface magnetization in the two- 
and the three-dimensional semi-infinite Ising models with 
rrii — 0.04. OT: ordinary transition, SP: special transition 
point. The grey line for the ordinary transition in three di- 
mensions is a fit to the predicted stretched exponential behav- 
ior JHJ at early times. The dashed lines indicate the power-law 
behavior taking place in the domain growth regime. 



two different regimes are observed. For x\ < Xi the mag- 
netization increases, whereas for x\ > xi it decreases, the 
decrease being nonalgebraic at early times. The pure sys- 
tem X\ — 1/2 is very close to the borderline value 0.53 
separating the two cases. This may explain the shal- 
low minimum observed after the very first time steps be- 
fore the magnetization increases again. In all cases a 
power-law behavior is encountered at later times, with 
an exponent that nicely agrees with the theoretical ex- 
pectation, see Table |H] The nonalgebraic decay in the 
cluster dissolution regime can in all cases be fitted by a 
stretched exponential. The values of the exponent k ex- 
tracted from these fits are listed in Table UTTI As for the 
three-dimensional Ising model at the ordinary transition 
we find a good agreement with the values obtained from 
the scaling theory. 

Finally, Figure is devoted to the two-dimensional 
Ising model with a ladder defect. As for all defect 
strengths J; one has xi < 1/2, the domain growth regime 
is the only one accessible within this model. Indeed, 
as shown in the Figure, we always observe an increas- 
ing magnetization. Again, the exponent extracted from 
these data perfectly agrees with the scaling prediction 
9i = (xi — xi)/z. It has to be noted that the same agree- 
ment is obtained for a chain defect when comparing the 
theoretical prediction with the numerically determined 
values given in Ref. |2!J 

Before turning to the autocorrelation functions, let us 
briefly discuss the robustness of the exponents 9i and k 
listed in Tables [n] and IIIII For this we investigated the 
ordinary transition in the pure semi-infinite systems by 
changing the strength of the surface couplings, J s . One 



10" 




FIG. 2: Relaxation of the surface magnetization in the 
Hilhorst-van Leeuwen model for different values of the surface 
scaling dimension x\. The value of the initial magnetization 
is vtii = 0.04. The grey lines for the cases x\ > Xi « 0.53 are 
obtained by fitting a stretched exponential to the early time 
data, see Eq. @. The expected power-law behavior is indi- 
cated by the dashed lines. The inset shows that for xi = 1 
ln(mi) is proportional to t K ^ z at early times, see Eq. 1A.1L 
as expected for the cluster dissolution regime. 



0.16 



0.10 - 




0.04 



FIG. 3: Relaxation of the surface magnetization in the Bariev 
model for different strengths of the defect coupling J;. The 
dashed lines indicate the expected power-law increase. 



then observes that the value of 9i is very robust against 
modifications of the surface coupling strength. In all 
studied cases a power-law behavior with a constant 9i 
sets in when entering the domain growth regime. When 
dealing with the three-dimensional system, the cluster 
dissolution regime is encountered at early times, yielding 
a stretched exponential behavior, with a crossover to a 
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power-law at later times. It is this crossover time which 
is strongly affected by the strength of the surface cou- 
plings. For weak surface couplings, as for example for 
Js = Jb/2, dissolution of clusters is extremely effective in 
the surface region, yielding a stretched exponential de- 
cay only at the very first time steps, the system crossing 
over to the DG regime very rapidly. This fact makes the 
extraction of k in this case very tedious. On the other 
hand, when increasing the strength of the surface interac- 
tions, one approaches the special transition point located 
at J s w 1.5J;,. As a result complicated crossover effects 
are encountered which again make a reliable determina- 
tion of k difficult. The value of k reported in Table ITTT1 
has been obtained for J s — Jt, where the discussed ef- 
fects are only of minor importance. Similar remarks also 
apply to the autocorrelation functions discussed in the 
following. 

C. Single spin autocorrelation function 

Following Eq. I|12|l the CD process should lead to an 
unusual stationary behavior of the single spin autocorre- 
lation at early times where t — s -c s. In our simulations 
we calculated the local quantity 

CM = | E ta(*)"i(«)> ( 2 6) 

i^surface 

with t > s, whereas S = L^ 1 is the size of the sur- 
face or interface and L is the linear extent of the system. 
The average is done over at least 10000 different realiza- 
tions of the thermal noise. The sum in Ij26|l is over all 
spins belonging to the surface or to the defect plane. In 
the present study we restrict ourselves to the short-time 
regime with t — s < 300. The dynamical scaling regime 
i, s, J - s > 1 has been analyzed in Ref. Izi . 

Our results for C(t, s) are summarized in Figures^and 
Eland in TablcHTll In FigureEJwe plot C(t, s) versus t — s 
for the pure semi-infinite cases, whereas the same is done 
in Figure in the Hilhorst-van Leeuwen model for two 
different values of the scaling dimension x\. In case of 
a stationary autocorrelation at short times the different 
curves corresponding to different waiting times should be 
indistinguishable. This is exactly what is observed for the 
cases with a CD regime where x\ > atr,. When x\ < Xi, 
however, the single spin autocorrelation function is not 
stationary but depends in a more complicated way on 
both times t and s and not merely on the time difference 
t — s. In the CD regime we again observe a stretched 
exponential decay as predicted, see Eq. i|12fl . This is 
illustrated in the inset of Figure 0Ji where we plot for the 
ordinary transition in three dimensions ln(C(t, s)) as a 
function of (t — s) K / z . Fitting a stretched exponential 
function to the short-time data we can extract the decay 
exponent k. The resulting values are listed in Table IIIII 
together with the estimates obtained from the relaxation 
measurements. A good agreement with the theoretical 



prediction is observed. 




FIG. 4: Single spin autocorrelation vs the time difference t — 
s for different waiting times s as obtained for the different 
pure semi-infinite systems: (a) ordinary transition and (b) 
special transition point in three dimensions, and (c) ordinary 
transition in two dimensions. When x\ > Xi a stationary 
behavior is observed with a stretched exponential decay at 
early times, see inset in (a). When xi < Xi the autocorrelation 
depends in a more complicated way on both times t and s. 




t-S t-S 

FIG. 5: Single spin autocorrelation as function of the time dif- 
ference t — s as obtained for the Hilhorst-van Leeuwen model 
with different waiting times. A stationary behavior is ob- 
served in case (b) where x\ > Xi. The grey line is obtained 
by fitting a stretched exponential to the short time data. 



D. Manifold autocorrelations 

The phenomenological theory presented in Section II 
yields very remarkable predictions regarding the surface 
manifold autocorrelation function in an inhomogeneous 
critical system. For a surface or a defect plane in two di- 
mensions the surface manifold M. has the dimensionality 
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d! = 1, whereas d! — 2 for a three-dimensional system. 
The surface manifold autocorrelation function is defined 
by 



G,(M) = (S(t)S(s)) 



(27) 



with S(t) = J2ieM the sum extending over all spins 

belonging to the manifold. 

For t — s scaling theory predicts a power-law increase 
of Gi(t,t) with an exponent {d! — xi — x)/z whenever 
d! > xi+x. When dl < xi+x, however, the surface man- 
ifold autocorrelation function should rapidly tend to a 
constant. Recalling the possible values of Xi and x in the 
models under investigation, we see that for the pure semi- 
infinite models and for the Bariev model Gi(t,t) should 
display a power-law increase in all cases. Only in the 
Hilhorst-van Leeuwen model should one enter the regime 
with d' < xi+x, with the amplitude value A = —0.75 (i.e. 
x,\ = 7/8) separating the two regimes. This limiting case 
x\ = 7/8 is shown in Figure^ by the grey line. Indeed, a 
completely different behavior is observed for smaller val- 
ues of x\ (lines above the grey line) when compared to 
larger values of x\ (lines below the grey line): in the for- 
mer case G\(t, t) increases with time whereas in the latter 
case Gi(t,t) rapidly saturates. In Figure 03 we compare 
the data with d' > xi + x to the expected power-law in- 
crease. For all studied cases we find excellent agreement 
with scaling theory. 

As noted in Section II scaling theory predicts for t > s 
in the DG regime with d' > xi + x a change in the 
form of the short distance expansion of the surface man- 
ifold autocorrelation function as compared with that for 
t = s. Looking at Eqs. (|18|) and (|19l) one remarks that 
the scaling combination s( x '~ x " z Gi(t, s)/G\(s, s) should 
then only be a function of t/s. As shown in Fig. 3 of Ref. 
Il8l plotting this scaling combination as a function of t/s 
yields a perfect data collapse for different values of the 
local scaling dimension xi and for different waiting times, 
thus demonstrating the expected t/s dependence. 

On the other hand in the CD regime with d! < xi + x 
the autocorrelation G\(t, s) should just be a function of 
the time difference t — s, see Eq. (|20|) . This stationary 
behavior is indeed observed, as shown in Figure \7\ for 
the Hilhorst-van Leeuwen model with x\ — 9/8. In the 
inset we verify the stretched exponential behavior 120|) 
by plotting ln(Gi) as a function of it — s) K / z for x\ = 1 
and x\ = 5/4. The values of n' obtained by fitting the 
short time data with this kind of function are gathered in 
Table llVl and compared with the theoretical expectation 
k' = (xi — x)d/(d— 1). For x\ = 5/4 we observe a simple 
exponential, in agreement with n' / z > 1. 



E. Persistence of surface manifold 

Finally, let us discuss the persistence of surface man- 
ifolds. The probability P pr {t) that the manifold magne- 
tization has not changed sign up to time t is of course 
directly related to the functional form of the surface 




0.006 



0.004 



100 



FIG. 6: (a) Surface manifold autocorrelation function for the 
Hilhorst-van Leeuwen model for various values of the scaling 
dimension: x\ — i/8 with i = 1, • • • , 10 (from top to bottom). 
The grey line is obtained for xi = 7/8 and separates the 
two different regimes with diverging Gi(t,t) for x\ < 7/8 
and saturating G\(t,t) for xi > 7/8. (b) Comparison of the 
increase of G\(t,t) for x\ < 7/8 with the expected power-law 
(dashes lines) in a log-log plot. From top to bottom: xi — i/8 
with i — 1, ■ • ■ , 6. 



manifold autocorrelation. Results for the CD regime are 
shown in Fig. 4 of Ref. Il8l Plotting ln(P pr ) as function 
of t K l z we observe straight lines that indicate that the 
probability is indeed governed by a stretched exponential 
decay whenever d! < xi+x. Values obtained for k! nicely 
agree with those obtained from the surface manifold au- 
tocorrelations, see Table llVl Figure Hand Table E| sum- 
marize our results for the DG regime with d' > xi + x. 
The same behavior is observed in the different models: 
the probability P pr {t) decays as a power of time. The 
values of the power-law exponent Q' pr extracted from the 
numerical data are gathered in Tabled No theoretical 
expression for Q' pr is available. 
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FIG. 7: Early time behavior of the surface manifold auto- 
correlation function in the CD regime as obtained for the 
Hilhorst-van Leeuwen model with xi = 9/8. In the inset the 
stretched (pure) exponential behavior is verified for two dif- 
ferent values of xi. 



TABLE IV: Numerically determined values of the exponent 
k! governing the stretched exponential behavior of the surface 
manifold autocorrelation function in the CD regime when d! < 
xi + x. The values have been extracted from (I) the behavior 
of surface manifold autocorrelations and (II) from the decay 
of the persistence probability P pr . The theoretical prediction 
is given by k' = (xi — x)d/(d — 1). 



Hilhorst-van Leeuwen model 



A 


Xl 


I 


II 


theoretical 


-0.75 


7/8 


1.41(2) 


1.41(2) 


1.41 


-1 


1 


1.63(2) 


1.61(2) 


1.65 


-1.25 


9/8 


1.86(2) 


1.84(2) 


1.88 


-1.50 


5/4 


2.03(2) 


2.03(2) 


2.04 



IV. DISCUSSION 

Nonequilibrium critical dynamics involves a time hori- 
zon, where the quench to the critical temperature is 
made, and this has a long-time, scale-free effect on the 
dynamical processes in the system. New nonequilibrium 
exponents enter into the theory: the relaxation of the 
magnetization and the spin-spin (and manifold) auto- 
correlations involve the scaling dimension of the initial 
magnetization, Xi , whereas persistence of different man- 
ifolds involve a new exponent pr . In a spatially inho- 
mogeneous system, such as at a free surface or at a de- 
fect plane, the local critical behavior is generally different 
from that in the bulk and one has to introduce local crit- 
ical exponents, such as xi for the local magnetization. In 
this paper we studied nonequilibrium critical phenomena 
in inhomogeneous systems, in which the dynamical be- 
havior is the result of an interplay between a spatial and 
a temporal inhomogeneity. Our main result is that in 



TABLE V: Numerically determined values of the manifold 
persistence exponent O' governing the power-law decay of the 
persistence probability P pr (t) in the domain growth regime. 
OT: ordinary transition, SP: special transition point. 



pure semi-infinite models 




Xl 


e' 


d= 2 


1/2 


3.0(2) 


d = 3, OT 


1.26 


4.6(3) 


d = 3, SP 


0.376 


0.61(2) 


Hilhorst-van Leeuwen model 


A 


Xl 


9' 


0.75 


1/8 


0.97(2) 


0.50 


1/4 


1.63(3) 


0.25 


3/8 


2.0(2) 


-0.25 


5/8 


3.9(3) 


-0.50 


3/4 


4.5(3) 


Bariev model 


Ji 


Xl 


6' 


0.2 


0.376 


2.3(1) 


0.4 


0.278 


1.75(5) 


0.6 


0.208 


1.27(3) 


0.8 


0.159 


0.90(2) 


1.0 


0.125 


0.71(2) 


1.4 


0.085 


0.55(1) 


1.8 


0.064 


0.40(2) 



such a system in early times two types of nonequilibrium 
processes can take place. Conventional domain growth 
phenomena is observed if the temporal surface is more 
disordered, than the spatial one. This is the rule in the 
bulk of the system. There is, however, a second, up to 
now unnoticed process: when the spatial surface is more 
disordered then the temporal one, then in early times lo- 
cal order is reduced and cluster dissolution takes place. 
This process is manifested by fast, stretched exponen- 
tial relaxation and by stationary, stretched exponential 
autocorrelations, which, however, involve a universal ex- 
ponent. We have shown by scaling theory and checked 
by numerical calculations that as far as relaxation and 
autocorrelations are concerned the local nonequilibrium 
processes can be fully characterized by the existing ex- 
ponents, xi, Xi and by the dynamical exponent, z, both 
in the DG and in the CD regimes. 

Our results can be applied or generalized to other 
types of systems. We mention that the cluster disso- 
lution scenario has recently been noticed in nonequilib- 
rium relaxation of reaction-diffusion systems, such as the 
contact process^ For another type of inhomogeneous 
systemsi2*i4 we mention inhomogeneities of geometrical 
origin, such as wedges, corners and cones. For example 
the local exponent at a corner, x c , depends on the open- 
ing angle ip, and in two dimensions it is related to 
the surface exponent as x c = x\-k /<p through conformal 
invariance^ Generally local order at a corner is weaker 
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FIG. 8: Manifold persistence probability P pr (t) as a function 
of t in the domain growth regime for various cases, a: defect 
line persistence in the Bariev model with J; = 1.8, b: sur- 
face persistence in the three-dimensional semi-infinite model 
at the special transition point, c: surface persistence in the 
Hilhorst-van Leeuwen model with A = 0.75, d: surface per- 
sistence in the Hilhorst-van Leeuwen model with A = 0.50, e: 
defect line persistence in the Bariev model with Ji = 0.4, f: 
surface persistence in the pure semi-infinite two-dimensional 
Ising model, g: surface persistence in the three-dimensional 
semi-infinite model at the ordinary transition. 

than at a free surface, therefore one expects that the clus- 
ter dissolution scenario often applies for such systems. 
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APPENDIX: NONEQUILIBRIUM DYNAMICS 
WITH CLUSTER DISSOLUTION 

Nonequilibrium dynamics of inhomogcneous sys- 
tems, in which the local (surface) scaling dimension, 
xi, is larger than the scaling dimension of the initial 
magnetization, Xi, is governed by cluster dissolution 
(CD)ii^ (Similar scenario holds for the surface manifold 
autocorrelation function for d' < xi + x.) In this case 
the nonequilibrium growth process, which goes in time 
as ~ t Xi ' z , is much weaker than critical local relaxation, 
which has a time dependence of t~ x '^ z . As a result 
typically no ordered domains are created in the surface 
region and the dynamics is governed by large correlated 
clusters which are present with a very small probability 
in the initial state. 



1. Relaxation 

Let us start with an initial state having a magnetiza- 
tion, rrii. The probability of having a cluster of linear 
size, I, is given by P(l) ~ m\ , where d is the dimension 
of the system. During the relaxation process the mass 
of a cluster is diluted by a factor of t( Xi ~ x, '' z , and the 
cluster is dissolved, if a domain wall of size V 1 " 1 is cre- 
ated in it. From the relation t(l) ( ' Xi ~ xl ^'l d ~ 1 = 0(1) we 
obtain for the typical time-scale of a cluster of linear size, 
I: t(l) ~ After time t > t(l) only those 

clusters exist, which have an original size larger than I. 
The magnetization, m(t), is just given by the contribu- 
tion of these large clusters: m(t) ~ J2i>i(t) -P(0^ ■ From 
this we obtain in leading order: 

\n(m(t)/ mi ) ~ -f/*, k= {X ^_ X i] d ■ (A- 1 ) 

Note, however, that the relaxation can not be faster, than 
in the non-critical case, when it is pure exponential, i.e. 
the exponent in Eci. (|A.lj) reads as: min(K/z,l). This 
result is compatible with the scaling form in Eq.© and 
0- 



2. Single spin autocorrelations 

The reasoning is similar as for relaxation and we repeat 
that for a completely uncorrelated initial state with zero 
magnetization the probability of having a large cluster of 
linear size, I, is given by P(l) ~ exp(— Al d ) and the cor- 
responding relaxation time is t r — t(l) ~ l( d - 1 ) z K x i~ x i) . 
The autocorrelation function is obtained by performing 
an average over the clusters: 

C(t) ~ J P(t r )exp(-T/t r )dt r ~ exp(-0W z ) , (A.2) 

in terms of r — t — s < s. If it happens that k/z > 
1, then the decay is pure exponential. Note that the 
autocorrelation function is stationary for early times and 
is compatible with Eq. l|12|l . 



3. Manifold autocorrelations 

For surface manifold autocorrelations the regime of 
cluster dissolution is valid for d' < Xi + x. In this 
case the dilution of ordered clusters goes in time as 
^{x-xi)/z^ As a consequence results for the single spin 
autocorrelation function can be easily translated, just in 
Ea. HA.2|) the decay exponent k has to be replaced by 
k' = (xi — x)d/ [d — 1). The corresponding results are in 
Eq.J2J. 
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